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Abstract 

We present a semi-analytic calculation of the tau-lepton flux emerging from the Earth, induced 
by the incident high energy neutrinos interacting inside the Earth for 10 5 < E u /GeV < 10 10 . We 
obtain results for the energy dependence of the tau-lepton flux coming from the Earth-skimming 
neutrinos, because of the neutrino-nucleon charged-current scattering as well as the resonant v e e~ 
scattering. We illustrate our results for several anticipated high energy astrophysical neutrino 
sources such as the AGNs, the GRBs, and the GZK neutrino fluxes. The tau lepton fluxes resulting 
from rock-skimming and ocean-skimming neutrinos are compared. Such comparisons can render 
useful information for the spectral indices of incident neutrino fluxes. 
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I. INTRODUCTION 



The detection of high energy neutrinos {E v > 10 5 GeV) is crucial to identify the extreme 
energy sources in the Universe, and possibly to unveil the puzzle of cosmic rays with energy 
above the GZK cutoff [lj]. These proposed scientific aims are well beyond the scope of the 
conventional high energy gamma-ray astronomy. Because of the expected small flux of the 
high energy neutrinos, large scale detectors (> 1 km 2 ) seem to be needed to obtain the first 
evidence. 

There are two different strategies to detect the footprints of high energy neutrinos. The 

first strategy is implemented by installing detectors in a large volume of ice or water where 

most of the scatterings between the candidate neutrinos and nucleons occur essentially inside 

the detector, whereas the second strategy aims at detecting the air showers caused by the 

charged leptons produced by the neutrino-nucleon scatterings taking place inside the Earth 

or in the air, far away from the instrumented volume of the detector. The latter strategy 

thus include the possibility of detection of quasi horizontal incident neutrinos which are also 

referred to as the Earth-skimming neutrinos. These neutrinos are considered to interact 

below the horizon of an Earth based surface detector. 

n 

The second strategy is proposed only recently [2j. The Pierre Auger observatory group 
has simulated the anticipated detection of the air-showers from the decays of r leptons[3(. 
The tau air-shower event rates resulting from the Earth-skimming tau neutrinos for different 
high energy neutrino telescopes are given in[J|. A Monte-Carlo study of tau air-shower event 
rate was also reported not long ago [5]. We note that Ref.jiJ does not consider the tau-lepton 
energy distribution in the z/ r -nucleon scattering, and only the incident tau neutrinos with 
energies greater than 10 8 GeV are considered. For Ref. [5(, we note that only the sum of 
tau air-shower event rates arising from different directions is given. Hence some of the 
events may be due to tau-leptons/neutrinos traversing a large distance. As a result, it is not 
possible to identify the source of tau-neutrino flux even with the observation of tau-lepton 
induced air-shower. 

In this work, we shall focus on the high energy Earth-skimming neutrinos and shall calcu- 
late the energy spectrum of their induced tau-leptons, taking into account the inelasticity 
of neutrino-nucleon scatterings and the tau-lepton energy loss in detail. Our work differs 
from Ref. p] by our emphasis on the Earth-skimming neutrinos. We shall present our results 
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in the form of outgoing tau-lepton spectra for different distances inside the rock, instead of 
integrating the energy spectra. As will be demonstrated, such spectra are insensitive to the 
distances traversed by the Earth-skimming v T and r. They are essentially determined by the 
tau lepton range. Because of this characteristic feature, our results are useful for setting up 
simulations with specifically chosen air-shower content detection strategy, such as detection 
of the Cherenkov radiation or the air fluorescence. 

We start with our semi-analytic description in Section II. The transport equations gov- 
erning the evolutions of neutrino and tau-lepton fluxes will be derived. Using these, we then 
calculate the tau-lepton flux resulting from the resonant u e e~ — > W~ — > u t t~ scattering. 
In Section III, we summarize our main results, namely the tau-lepton energy spectra due 
to neutrino-nucleon scatterings. The implications of our results will be discussed here also. 
In particular, we shall point out that the ratio of tau-lepton flux induced by rock-skimming 
neutrinos to that induced by ocean-skimming neutrinos is sensitive to the spectral index of 
the incident tau-neutrino flux. In Section IV, we discuss some prospects for possible future 
observations of the associated radiation from these tau leptons. 



II. TAU LEPTON ENERGY SPECTRUM 



Let us begin with the transport equations for tau neutrinos and tau leptons. Considering 
only the neutrino-nucleon scatterings, we have 

dF Vr (E,X) F Ut (E,X) , A dy _ _d<r; 



dX ~ \ Ut (E) 



■J^ fVlnax Ay Apt 

Jy l . 1 ~ y a y 



and 



dF T (E, X) _ F T (E, X) F T (E, X) d [(a(E) + (3(E)E)F T (E, X)} 



OX X C T C (E) pd T (E) dE 

T max &(j VtN ^ t y , , . 

+ Un / ii:^^^) — — ^' ^' ® 

where the number of target nucleons per unit medium mass, and p is the mass density 
of the medium. The a*' 2 ' 3 are defined as a(u T + N — > v T + Y), T(t — > u T + Y)/cpriN, and 
a{r + N — > u T + Y) respectively. The quantity X represents the slant depth traversed by 
the particles, i.e., the amount of medium per unit area traversed by the particle (and thus 
in units of g/cm 2 ). The \ u , d T , and A^ c represent the v T interaction thickness, the tau- 
lepton decay length, and the tau-lepton charged- current interaction thickness respectively, 



with, say, A" 1 = unct u n and d T = CT T E/m T . The E y is equal to E/(l — y), where y is 
the inelasticity of neutrino-nucleon scatterings, such that the initial and final-state particle 
energies in the differential cross sections da l u (y, E y )/dy and da VT N^ T y(y, E y ) / dy are E / (1—y) 
and E respectively. The limits for y, y l min and y^^ depend on the kinematics of each process. 
Finally, the energy-loss coefficients a(E) and (3(E) are defined by — dE/dX = a(E)+(3(E)E 
with E the energy being the tau lepton. An equation similar to Eq. (J2J) in the context of 
atmospheric muons was found in Ref. Q]. 

As mentioned before, Eqs. (0) and (J2J) only take into account neutrino-nucleon scatterings. 
It is of interest to calculate the tau-lepton fluxes produced by the Glashow resonance 0, Ea], 
namely via v e e~ — > W —> v t t~ also. The transport equation for v e then reads: 

a^JE.X) _ F,JE,X )+nN r -^ F ^ x) ^r^ Eyl (3) 



9X K.(E) "J ymm l-y " ' dy 

Similarly, the corresponding equation for the tau-lepton flux is given by 

dF T (E,X) F T (E,X) F T (E,X) 
OX \?(E) pd T (E) 

ft/max A.. A (J _ _ _ _ 

+ n e ^4(^,1) ^ (y,E y ), (4) 

where n e is the number of target electrons per unit medium mass. 

Before solving the above coupled transport equations, it is essential to know the energy- 
loss coefficients a(E) and (3(E). As pointed out before |9(, the coefficient a(E) is due to the 
energy loss by the ionization 10( , while (3(E) is contributed by the bremsstrahlung 1 JJ , the 



:d Dy t 

0,Q. 



e + e~ pair production |12j and the photo-nuclear processes |9l. It is understood that the 
contribution by a(E) becomes unimportant for E > 10 5 GeV. The coefficient (3(E) can be 
parameterized as (3(E) = (1.6 + 6(i?/10 9 GeV) ' 2 ) x 10~ 7 g -1 cm 2 in the standard rock for 
10 5 < E/GeV < 10 12 . 

It is of interest to check the tau-lepton range given by our semi-analytic approach. To 
do this, we rewrite Eq. (J2j) by dropping the neutrino term, i.e., 

dF T (E,X) F T (E,X) F T (E,X) d[^(E)F T (E,X)] 

dX \?(E) pd T (E) dE U 

with 'j(E) = a(E) + (3(E)E. One can easily solve it for F T (E, X): 

■x 



F T (E,X) = F T (E,0)exp 



I dT ^ E) -^)->m) 



(6) 
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where E = E(X;E) with dE/dX = 7(E) and E(0; E) = E. To calculate the tau-lepton 
range, we substitute F T (E,0) = 5(E — E ). The survival probability P(E ,X) for a tau- 
lepton with an initial energy E at X = is 



P(B„,A-) = ^exp 
J{E ) 



r dr ( i{E ) — k — L-^i 

Jo \ pd T (E ) \?(Eo)J 



(7) 

L^o 

where E = E(X;E ) with dE /dX = —7(^0) and E (0; E ) = E . The tau-lepton range 
is simply 

POO 

R T (E )= / dXP(E ,X). (8) 
Jo 

For £ = 10 9 GeV, we find that R T = 10.8 km in the standard rock (Z = 11, A = 22) 
while R T = 5.0 km in the iron. Both values are in good agreement with those obtained 



by Monte- Carlo ca 



culations[9]. To compare the tau-lepton ranges, we have followed the 



by requiring the final tau-lepton energy E(X; E ) to be greater than 



convention in Ref. 
50 GeV. 

It is to be noted that we obtain R T by using the continuous tau-lepton energy- loss 
approach, rather than stochastic approach adopted in Ref.j^|. In the muon case, the con- 
tinuous approach to the muon energy loss is known to overestimate the muon rangejl^. 
Such an overestimate is not significant in the tau-lepton case, because of the decay term in 
Eq. 0. In fact, tau lepton decay term dictates the tau range in the rock until E T > 10 7 
GeV. Even for E T > 10 7 GeV, the tau lepton range is still not entirely determined by the 
tau-lepton energy loss. Hence different treatments on the tau-lepton energy loss do not lead 
to large differences in the tau-lepton range, in contrast to the case for the muon range. Our 
result for the tau-lepton range up to 10 12 GeV are plotted in Fig. ^ This is an extension 
of the result in Ref.Jsj], where the tau-lepton range is calculated only up to 10 9 GeV. Our 
extension is seen explicitly in the addition of charged- current scattering term on the R.H.S. 
of Eq. (jSJ). This term is necessary because 1/A^ C becomes comparable to l/pd T in the rock 
for E > 10 10 GeV; whereas one does not need to include the contribution by the tau-lepton 
neutral-current scattering, since such a contribution can not compete with the last term in 
Eq. (jHJ) until E > 10 16 GeV[3|. We remark that our extended results for R T are subject to 
the uncertainties of neutrino-nucleon scattering cross section at high energies. We use the 
CTEQ6 parton distribution functions ^lj]|] in this work, and at the high energy (the small 
x region, namely for x < 10~ 6 ), we fit these parton distribution functions into the form 
proportional to £~ L3 as a guide. 



Having checked the tau-lepton range, we now proceed to calculate the tau-lepton flux. It 
is instructive to begin with the simple case: the v e e~ resonant scattering. It is well known 
that HQ 

^-.r.^-)^. ^^^ , (9) 

with s = 2m e Ev e and 1/a-da/dz = 3(1 — z) 2 , where z = E T /E Pe . We shall only focus on those 
zVs for which E Pe satisfies the resonance condition, i.e., E Pe fa E R = m 2 v /2m e . It is clear 
from Eq. (jlj) that F T (E, X) only depends on F 9e (E R , X), because of the narrow peak nature 
of v e e~ scattering cross section. One also expects that F T (E,X) is only significant for E 
around the resonance energy E R . In this energy region, one may neglect the first term on the 
R.H.S. of Eq. in comparison with the second term. In the narrow width approximation, 
the last term in Eq. (J3| can be recast into — E / E R ) 2 (iiT w / L R m w )F De (E R , X), where T w 
is the width of the W boson while L R is the interaction thickness of the resonant v e e~ — > W~ 
scattering (see Appendix A for details). The tau lepton flux can be readily obtained once 
Fp e (E R , X) is given. We observe that the regeneration term in Eq. © (second term on the 
R.H.S.) can be neglected as it is necessarily off the W boson peak. Hence, we easily obtain 
F p (E R ,X) = exp(—X/L R )Fp (E R ,0). Substituting this expression into Eq. (jU), we obtain 

IS— (IH- ->(-£)■ ™ 

in the limit X 3> pd T . The pre-factor 3.3 -10 4 is obtained by assuming a standard-rock 
medium. In water it becomes 1.4 • 10~ 4 . It is to be noted that E < E R in the above 
equation. We shall see later that the contribution to F T (E,X) by the VT-resonance is 
negligible compared to that by the v T — N scattering. 

Let us now turn to the case of tau-lepton production by the v T — N charged-current 
scattering. The tau-lepton flux can be calculated from Eqs. (0) and J21 once the incoming 
v T flux is given. The v T flux can be obtained by the following ansatz|l6j]: 

F„( £[ A>F„( f ,0) np (-^), (11) 

where A V (E,X) = X V (E)/ (1 — Z V (E,X)), with the factor Z U (E,X) arising from the regen- 
eration effect of the v T flux. On the other hand, the tau-lepton flux is given by 

F T (E,X) = [ dTG v (E, T) 



x exp 



(12) 
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with E = E(X - T; E), and 

G v (E,X)=n N -^-E u (E y ,X) aa ^ TY (y,E y ). (13) 

It is easy to see that the factor Z l/ (E, X) enters into the expression for F T (E, X) through the 
function G„(E, X). Similarly, Z V (E, X) also depends on F T (E, X). It is possible to solve for 
Z U (E,X) and F r (E,X) simultaneously by the iteration method[16]. The details are given 
in Appendix B. 

III. RESULTS AND DISCUSSION 

In the following, we show the tau-lepton fluxes resulting from three kinds of diffuse 
, tr oph y siealne» t lo fluxes: the AGN0, GRB0 and GZK 191 neutrino fluxes. In these 
representative models, F Vt arises because of neutrino flavor mixing The p7 interactions 
are the source of intrinsic F v , and F Vt = 1/2 • F v because of (two) neutrino flavor oscillations 
during propagation. Our convention for F Ut is that F Ut = dN Ur / d(log w E) in the unit of 
cm -2 s _1 sr _1 . The same convention is used for the outgoing tau lepton fluxes. 

In Fig. |2l we show the outgoing tau-lepton energy spectra resulting from the propagation 
of incident AGN neutrinos inside the rock (p = 2.65 g/cm 3 ) for X/p = 10 km, 100 km 
and 500 km respectively. It is interesting to see that the tau-lepton energy spectra remain 
almost unchanged for the above three different slant depths/matter density ratio values. This 
feature can be understood by two simple facts. First of all, the neutrino-nucleon charged- 
current interaction length, which is related to the interaction thickness by Acc = P^cc? is 
given by Zee = 2- 10 4 km ( - g ^ cm )( 10 <^ eV )~°' 363 - Secondly, the tau leptons, which eventually 
exit the Earth, ought to be produced within a tau-lepton-range distance to the exit point. 
For a tau-lepton produced far away from the exit point, it loses energy and decays before 
reaching to the exit point. Hence the tau-lepton flux is primarily determined by the ratio of 
tau-lepton range to the charged current neutrino-nucleon interaction length. The total slant 
depth X which the tau-neutrino (tau-lepton) traverses inside the Earth is then unimportant, 
unless X is large enough such that the tau neutrino flux attenuates significantly before 
tau-neutrino is converted into the tau-lepton. We note that the typical energy for the 
AGN neutrinos, in which this flux peaks, is between 10 5 and 10 8 GeV. The corresponding 
neutrino-nucleon neutral current interaction length then ranges from 42, 000 km down to 
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3,400 km, given /nc = 2.35 • Zee- Hence, even for X/p as large as 500 km, the attenuation 
of the tau neutrino flux is negligible. This explains the insensitivity of tau-lepton flux with 
respect to our chosen Xj p values for the AGN case. The situation is rather similar for the 
tau-lepton flux resulting from the GRB tau neutrinos (see Fig. On the other hand, a 
slight suppression is found for the GZK case at E T > 10 9 GeV as one increases X/p from 
10 km to 500 km (see Fig. HJ. This is because the typical GZK tau neutrino flux peaks for 
energy range between 10 7 and 10 10 GeV, which corresponds to attenuation lengths ranging 
from 7, 800 km down to 640 km. One notices that 640 km is rather close to the distance 500 
km which we choose for Xj p. Hence a slight suppression on the tau-lepton flux occurs for 
X/p = 500 km. 

We have compared our AGN-type tau-lepton flux with that obtained by Monte-Carlo 
simulations, adopting stochastic approach for the tau-lepton energy loss |21l| . The two tau- 
lepton fluxes agree within ~ 10%. This is expected since the tau-lepton range obtained 
by the above two approaches agree well, as pointed out before. It is easily seen from Fig. 
El to Fig. 13] that the AGN case has a largest tau-lepton flux between 10 6 and 10 8 GeV. 
Since the resonant v e — e~ scattering cross section peaks at E v = 6.3 • 10 6 GeV, it is of 
interest to compare the integrated tau-lepton flux resulting from this scattering to the one 
arising from neutrino-nucleon scattering. For the former case, we integrate the tau-lepton 
energy spectrum from 10 6 GeV to 6.3 • 10 6 GeV, and obtain $^ = 0.08 km~ 2 sr -1 yr _1 . 
For neutrino-nucleon scattering, we find that = 2.2 km -2 sr~ 1 yr _:l by integrating the 
corresponding tau-lepton energy spectrum from 10 6 GeV to 10 7 GeV. The detailed results for 
are summarized in Table HJ The entries in the Table entitled full are obtained using the 
F T obtained in this work, whereas the approximated values entitled approx are obtained by 
following the description given in Ref. ^], which uses a constant (3 and a constant inelasticity 
coefficient for u T N scattering. We remark that the authors in Ref. 4J have taken E to be 
greater than 10 8 GeV. Hence the integrated fluxes in the column approx with energies less 
than 10 8 GeV are taken as extrapolations. Thus, one should compare the two integrated 
fluxes only for E > 10 8 GeV. One can see that the two integrated fluxes seem to agree for 
E > 10 8 GeV. Besides the integrated fluxes for E > 10 8 GeV, we also obtain integrated 
tau-lepton fluxes for 10 6 < E/GeV < 10 8 . It is easily seen that, in this energy range, the 
integrated tau-lepton flux from Earth-skimming AGN neutrinos is relatively significant . 

It is possible that the tau-neutrino skims through a part of the ocean in addition to the 
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Earth before exiting the interaction region 22] . Hence, it is desirable to compare the resulting 
tau-lepton fluxes as the tau neutrinos skim through mediums with different densities, while 
the slant depths of mediums are held fixed as an example. As stated before, the tau- 
lepton flux is essentially determined by the probability of v T N charged-current interaction 
happening within a tau-lepton range. Furthermore, from Fig. ^ it is clear that the tau- 
lepton range equals to the tau-lepton decay length for E T less than 10 7 GeV. One therefore 
expects F™ ck (E,X)/F™ tei (E,X) = p rock / p water for E T < 10 7 GeV. This is clearly seen to be 
the case from Fig. Eland Fig. H as we compare F T rock with F r water (£, X) for X = 2.65 • 10 6 
g/cm 2 and X = 2.65 • 10 g/cm 2 respectively. For E T > 10 7 GeV, the tau-lepton range has 
additional dependencies on the mass density and the atomic number of the medium. Hence 
the ratio F r rock (£, X)/F™ ter (E, X) starts deviating from pP^/p™**. It is worthwhile to 
mention that the tau-lepton flux ratios for AGN and GRB cases behave rather similarly. On 
the other hand, the ratio in the GZK case has a clear peak between 10 7 5 < E/GeV < 10 8,5 . 
Such a peak is even more apparent for the slant depth X = 2.65 -10 7 g/cm 2 . The appearance 
of this peak has to do with the relatively flat behavior of the incident GZK neutrino spectrum, 
while the position of this peak is related to the energy dependencies of the tau-lepton range 
and the neutrino-nucleon scattering cross sections. We have confirmed our observations by 
computing the flux ratios with simple power-law incident tau-neutrino fluxes. The above 
peak in the tau-lepton flux ratio implies the suppression of tau-lepton events from ocean- 
skimming neutrinos compared to those from rock-skimming neutrinos. As stated earlier, 
the suppression of ocean-skimming neutrinos is related to the spectral index of the incident 
neutrino flux. It is therefore useful to perform a detailed simulation for it [^J. Such a detailed 
study is needed because the slant depths traversed by the above two kinds of neutrinos are 
generally different. 

IV. PROSPECTS FOR POSSIBLE FUTURE OBSERVATIONS 

To observe the above tau leptons, the acceptance of a detector must be of the order of 
~ km 2 sr. For AGN neutrinos, the tau-lepton energy spectrum peaks at around 10 7 to 10 8 
GeV, which is below the threshold of a fluorescence detector, such as the High Resolution 



Fly's Eye (HiRes) 2M- Also, these tau leptons come near horizontally. At present, it seems 



very difficult to construct a ground array in vertical direction. A Cherenkov telescope seems 
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to be a feasible solution. In this context, NuTel collaboration is developing Cherenkov 



22 1. However, because of the 



telescopes to detect the Earth-skimming high energy neutrinos 
small opening angle of Cherenkov light cone and only a 10% duty cycle (optical observations 
are limited to moonless and cloudless nights only), such a detector must cover very large 
area and field of view. A potential site for NuTel is at Hawaii Big Island, where two large 
volcanos, namely Mauna Loa and Mauna Kea, could be favorable candidates for high energy 
neutrinos to interact with. For a detector situated on top of Mount Hualalai and to look 
at both Mauna Kea and Mauna Loa, the required angular field of view is ~ 8° x 120°. 
Furthermore, this telescope should have an acceptance area larger than 2 km 2 sr so as to 
detect more than one event per year. 

Concerning the GZK neutrinos, we note that the recent observation of ultra high energy 
cosmic rays by HiRes seem to be consistent with the GZK cutoff. Therefore a future obser- 
vation of GZK tau neutrinos shall provide a firm support to GZK cutoff. In particular, the 
slight pile up of tau leptons between 10 8 GeV to 10 9 GeV, induced by the Earth-skimming 
high energy GZK neutrinos, should be a candidate signature for GZK neutrinos. The inte- 
grated tau-lepton flux in this energy range is approximately 0.08 km -2 sr -1 yr -1 . To detect 
one event per year from this flux, the acceptance of a detector must be larger than 120 
km 2 sr, for a fluorescence detector (assuming a duty cycle of 10% ). Although HiRes can 
reach 1000 km 2 sr at energy > 3 • 10 9 GeV, it would be a technical challenge to lower down 
the threshold to 10 s GeV. Using a system similar to HiRes, the Dual Imaging Cherenkov Ex- 
periment (DICE) was able to detect Cherenkov light from extensive air-showers at energy as 
low as 10 5 GeV[25]. However, the field of view of DICE is also quite small, and thus several 
Cherenkov telescopes would be needed. An alternative method is a hybrid detection of both 
Cherenkov and fluorescence photons That is, a detector similar to HiRes, which looks 
at both land and sea and detects both Cherenkov and fluorescence photons, may observe 
the associated signal of GZK neutrinos. 

In summary, we have given a semi-analytic treatment on the problem of simultaneous 
propagation of high energy tau neutrinos and tau leptons inside the Earth. Our treatment 
explicitly takes into account the inelasticity of neutrino-nucleon scatterings as well as 
the tau-lepton energy loss. We specifically considered the Earth-skimming situation and 
provided detailed results for the energy dependencies of emerging tau-lepton fluxes resulting 
from a few anticipated astrophysical neutrino fluxes. The effect of matter density on the 



10 



tau-lepton flux is also studied. Such an effect is found to be related to the spectrum index of 
incident neutrino flux. Our treatment thus provides a basis for a more complete and realistic 
assessment of high-energy- neutrino flux measurements in the under-construction/planning 
large neutrino telescopes. 

Acknowledgements 

We thank N. La Barbera for communicating to us his Monte- Carlo-based results. HA 
thanks Physics Division of NCTS for support. MAH is supported by Taiwan's Ministry of 
Education under Research Excellence Project on Cosmology and Particle Astrophysics: Sub- 
project //with the grant number 92-N-FA01-1-4-2. FFL, GLL, JJT and TWY are supported 
by National Science Council of Taiwan under the grant numbers NSC91-2112-M009-019 and 
NSC91-2112-M-001-024. 



[1] S. W. Barwick, Astroparticle physics with high energy neutrinos, lectures presented at 28th 
SLAC Summer Institute on Particle Physics: Neutrinos from the Lab, the Sun, and the Cosmos 
(SSI 2000), Stanford, California, USA, 14-25 Aug 2000; H. Athar, |a7Xiv:hep-ph/0209T30 For 
a review article, see H. Athar, arXiv:hep-ph/0212387 

[2] G. Domokos and S. Kovesi-Domokos, |arXiv:hep-ph/ 980 1362] |arXiv:hep-ph/980522T| See also, 



D. Fargion, Astrophys. J. 570, 909 (2002) ar Xiv:astro-ph /0002453 . 
[3] X. Bertou, P. Billoir, O. Deligny, C. Lachaud and A. Letessier-Selvon, Astropart. Phys. 17, 

183 (2002) |arXiv:astro-ph/0104452| . 
[4] J. L. Feng, P. Fisher, F. Wilczek and T. M. Yu, Phys. Rev. Lett. 88, 161102 (2002) 

|arXiv:hep-ph/01050"67] . 
[5] S. Bottai and S. Giurgola, Astropart. Phys. 18, 539 (2003) |arXiv:astro-ph/0205325| . 
[6] L. V. Volkova, G. T. Zatsepin and L. A. Kuzmichev, Sov. J. Nucl. Phys. 29, 645 (1979) [Yad. 

Fiz. 29, 1252 (1979)]. 
[7] S. L. Glashow, Phys. Rev. 118, 316 (1960). 

[8] For a recent discussion see, H. Athar and G. -L. Lin, Astropart. Phys. 19, 569 (2003) 
arXiv:hep-ph/0108204 , and references cited therein. 



11 



[9] S. I. Dutta, M. H. Reno, I. Sarcevic and D. Seckel, Phys. Rev. D 63, 094020 (2001) 
|arXiv:hep-ph/0012350] . 
[10] B. Rossi, High Energy Particles (Prentice Hall, Englewood Cliffs, NJ, USA, 1952). 
fill A. A. Petrukhin and V. V. Shestakov, Can. J. Phys. 46, S377 (1968). 
[12] R. P. Kokoulin and A. A. Petrukhin in Proceedings of the XII International Conference on 

Cosmic Rays (Hobart, Tasmania, Australia, 1971), Vol. 6. 
[13] L. B. Bezrukov and E. V. Bugaev, Sov. J. Nucl. Phys. 33, 635 (1981) [Yad. Fiz. 33, 1195 
(1981)]. For a recent discussion, see E. V. Bugaev and Y. V. Shlepin, Phys. Rev. D 67, 034027 



[14 
[15 

[16 
[17 



[18 

[19 

[20 

[21 
[22 



(2003) |arXiv:hep -ph/0203096 . 

P. Lipari and T. Stanev, Phys. Rev. D 44, 3543 (1991). 

J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207, 
012 (2002) |arXiv:hep-ph/0201195|. 



V. A. Naumov and L. Perrone, Astropart. Phys. 10, 239 (1999) arXiv:hep-ph/9804301 . 
A. Neronov, D. Semikoz, F. Aharonian and O. Kalashev, Phys. Rev. Lett. 89, 051101 (2002) 
|arXiv:astro- ph/0201410] ; O. E. Kalashev, V. A. Kuzmin, D. V. Semikoz and G. Sigl, Phys. 
Rev. D 66,063004 (2002) |arXiv:hep-ph/0205050| . 

E. Waxman and J. N. Bahcall, Phys. Rev. Lett. 78, 2292 (1997) |arXiv:astro-ph/970123l"1 ; 
ibid., Phys. Rev. D 59, 023002 (1999) |arXiv:hep-ph7 9807282 . 

V. S. Berezinsky and G. T. Zatsepin, Phys. Lett. 28B, 423 (1969). For a recent discussion, see 
R. Engel, D. Seckel and T. Stanev, Phys. Rev. D 64, 093010 (2001) |arXiv:astro-ph/0101216| . 



H. Athar, arXiv:hep-ph/0210244; H. Athar, K. Cheung, G. -L. Lin and J. -J. Tseng, Astropart. 
Phys. 18, 581 (2003) |arXiv:hep-ph/0112222| . 
N. La Barbera, private communications. 

G. W. Hou and M. A. Huang, Expected performance of a neutrino telescope for seeing 
AGN/GC behind a mountain, talk given in 1st NCTS Workshop on Astroparticle Physics, 
Taiwan, 6-9 Dec 2001, published in its proceedings, pp. 105-116, Eds. H. Athar, G. -L. 
Lin and K. -W. Ng, World Scientific, Singapore (2002) a rXiv:astro-p h/0204145 . See also 



http://hepl.phys.ntu.edu.tw/nutel/ 



[23] M. A. Huang, G. -L. Lin and J. -J. Tseng, work in progress. 

[24] T. Abu-Zayyad et al. [HiRes-MIA Collaboration], Astrophys. J. 557, 686 (2001) 



arXiv:astro-ph/0010652 . 

12 



[25] S. P. Swordy and D. B. Kieda, Astropart. Phys. 13, 137 (2000) |arXiv:astro- ph/990938l"1. 
[26] Z. Cao, M. A. Huang, P. Sokolsky, and R. W. Springer, Proc. of the 28th Int'l Cosmic Ray 
Conf., Tsukuba, Japan, July- August 2003, in preparation. 



13 



APPENDIX A: THE CONTRIBUTION FROM RESONANT P e e SCATTERING 



The transport equations for v e and the tau lepton are given by Eqs. © and (HJ). For 
convenience, let us write 1 — y — z. The last term in Eq. (j3J) can be simplified using 

and the narrow-width approximation 

1 7Tl\yTw ^ ^2 

7T (s - m 2 w ) 2 + m 2 w T 2 w 

We arrive at 



S(s-m 2 w ). (A2) 



dF T (E,X) _ F T (E,X) 1 / E\ 2 (7rT w \ 

where Er = m 2 v /2m e is the v e energy such that the W boson is produced on-shell in the 
v e e~ scattering. The Lr = l/n e a Pee -^ w - is the interaction thickness for such a scatter- 
ing. To solve for F T (E,X), we need to input F Pe (Er, X). Obviously, the v e flux at the 
resonant-scattering energy Er is mainly attenuated by the resonant scattering itself. Hence 
F De (E R ,X) = exp(—X/L R )Fp e (ER,0). Substituting this result into Eq. (jA3|) . we obtain 

^-K i -s) , (^)^"- p (-^) 

i i 



X 



x / dZ exp 



Z 



(A4) 



o LK pd T (E) L Rj 

The integration over Z can be easily performed. In practice, it is obvious that X ^> pd T (E). 
In this limit, we have 

Let us consider standard rock as the medium for v e e~ scattering, we then have p/Lr = 
n e pa 9ee -^w-- Given a„ ee -_+w- = 4.8 • 10~ 31 cm 2 at the W boson mass peak, and n e p = 
2.65 x 6.0/2 x 10 23 /cm 3 in the standard rock, we obtain p/Lr = (26 km) -1 . Furthermore, 
we can write d T (E) = 49km x (i?/10 6 GeV). We then obtain the following ratio 



F& e (E R ,0) \E R J \ Er) \ L R 

This is the result given by Eq. (fTU|) in the main text. 
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APPENDIX B: THE ITERATION METHOD FOR OBTAINING THE Z V (E,X) 

AND THE F T (E,X) 



The evolution for F Vt is given by Eq. With the ansatz 

F Ur {E,X) = iU£,0)exp (- -^g-^ ) , (Bl) 
we obtain the following equation for Z V (E,X): 

XZ V {E,X) = [ X dX' f -^l^^-ex P [-X'D u (E,E y ,X')]^ T C (y,E) 
Jo Jo 1-V { F£>(E) 

FSHE) \pd T (E)) eKP W(E,X'))® Ay ^ > 
F T (E„X>) (\„(E)\ ( X' \ co 1 

where Fv~}(E) = F Ut (E,0), while $^ T C , and <J>^ are respectively given by 

da UT T^v T Y 



^ = ^ ^F7^ > (B3) 



^ = 5^ ' (B4) 

^ = f^^("^ (B5) 
with ut the number of targets per unit mass of the medium, and 

For the simplicity in notations, we take the lower and upper limits for the y integration to be 
and 1 respectively. In reality, the limits depend on the actual kinematics of each process. 
One may impose these limits in the functions $^ t c , and 

To perform the iteration, we begin by setting Z v (q) = 0. In this approximation, we have 

F Ut{0) (E,X) = F UT (E,0)exp (- — |— ) . (B7) 
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Substituting F Ut ^(E, X) into Eq. (fT2j). we obtain the lowest order v T flux, F T ^(E, X). 
The first iteration for Z u , denoted by Z u ^ is calculable from Eq. (jB2|) by substituting 
F UT (p)(E, X), F T (fy(E,X), and Z u ^ into the R.H.S. of this equation. From Z u ^, we can 
then calculate F Ut ^(E, X) and F T n\(E,X), which corresponds to the results presented in 
this paper. We have checked the convergence of iteration procedure and have found negligible 
differences between Z u {n\ and Z v ^ and their associated v T and r fluxes. 

The value of Z v depends on the spectrum index of the neutrino flux, since it effectively 
gives the regeneration effect in the neutrino-nucleon scattering. In general, a flatter neutrino 
spectrum implies a larger Z v . The Z v is however not sensitive to the slant depth X. In the 
case of GRB neutrinos, where the flux decreases as E~ 2 for E v < 10 7 GeV, and decreases 
as E~ 3 for energies greater than that, we obtain Z GRB « 0.2. For the AGN neutrino, Z^ GN 
changes from 0.96 to 0.35 as E v runs from 10 5 GeV to 10 6 GeV. In this energy range, the 
neutrino flux decreases slower than E~ ^. For E u greater than 10 8 GeV, Z^ GN drops below 
0.2 as the neutrino flux spectrum begins a steep fall. The values for Z GZK also follow the 
similar pattern. 
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TABLE I: Comparison of the integrated tau-lepton flux (km _2 yr _1 sr _1 ) in different energy bins for 
the AGN, the GRB and the GZK neutrinos without and with approximation (see text for details). 
The distance traversed is taken to be 10 km in rock here. For 10 9 < E/GeV < 10 10 , the incident 
AGN neutrino flux is too small so that its induced tau-lepton flux is not shown. 



Energy Interval 


AGN 


GRB 


GZK 


full 


approx 


full 


approx 


full 


approx 


10 6 < E/GeV < 10 7 


2.23 


2.12 


9.63 • 10~ 3 


1.05 • 10~ 2 


7.38 • 10~ 5 


2.08 • 10" 5 


10 7 < E/GeV < 10 8 


4.89 


5.12 


7.12 • 10~ 3 


6.82 • 10" 3 


1.14 • lO" 2 


1.90 • 10^2 


10 s < E/GeV < 10 9 


1.95 • 10- 1 


1.52 • 10" 1 


5.39 • 10~ 4 


4.63 • 10" 4 


8.17-10" 2 


8.47 • 10^2 


10 9 < E/GeV < 10 10 






1.13 • 10" 5 


1.24 • 10~ 5 


3.31 • 10^2 


3.52 • 10-2 



17 



1 I ' 1 1 1 1 r 



10 2 N" 




10" 



Tau lepton decay length 

Tau lepton range in water 

- - - Tau lepton range in rock 



7 8 9 

log(E/GeV) 



10 11 



12 



FIG. 1: The tau- lepton range in rock and in water using Eq.[Hj and the tau- lepton decay length 
d T in km as a function of tau-lepton energy in GeV. 
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log(E / GeV) 

FIG. 2: The tau-lepton energy spectrum induced by the AGN neutrinos in rock for three different 
X/p ratio values (see text for more details). The incident tau-neutrino flux is shown by the thin 
solid line. 
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log(E / GeV) 

FIG. 3: The tau-lepton energy spectrum induced by the GRB neutrinos in rock for three different 
X/p ratio values (see text for more details). The incident tau-neutrino flux is shown by the thin 
solid line. 
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log(E / GeV) 

FIG. 4: The tau-lepton energy spectrum induced by the GZK neutrinos in rock for three different 
X/p ratio values (see text for more details). The incident tau-neutrino flux is shown by the thin 
solid line. 
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log(E / GeV) 

FIG. 5: The ratio of F T in rock and water induced by the AGN, the GRB and the GZK neutrinos 
for X = 2.65 • 10 6 g/cm 2 . 
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FIG. 6: The ratio of F T in rock and water induced by the AGN, the GRB and the GZK neutrinos 
for X = 2.65 • 10 7 g/cm 2 . 
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